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ABSTRACT 

We present mid-infrared observations of active galactic nnclei (AGN) in the 
Great Observatories Origins Deep Snrvey (GOODS) fields, performed with the 
Spitzer Space Telescope. These are the deepest infrared and X-ray fields to date 
and cover a total area of ~0.1 sqnare degrees. AGN are selected on the basis 
of their hard (2-8 keV) X-ray emission. The median AGN infrared Inminosity 
is at least 10 times larger than the median for normal galaxies with the same 
redshift distribntion, snggesting that the infrared emission is dominated by the 
central nnclens. The X-ray to infrared Inminosity ratios of GOODS AGN, most 
of which are at 0.5<z<1.5, are similar to the valnes obtained for AGN in the 
local Universe. The observed infrared flnx distribntion has an integral slope 
of ~1.5 and there are 1000 sonrces per sqnare degree brighter than ~50 /iJy 
at ~3-6 microns. The connts approximately match the predictions of models 
based on AGN nnification, in which the majority of AGN are obscured. This 
agreement conhrms that the faintest X-ray sources, which are dominated by the 
host galaxy light in the optical, are obscured AGN. Using these Spitzer data, the 
AGN contribution to the extragalactic infrared background light is calculated 
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by correlating the X-ray and infrared catalogues. This is likely to be a lower 
limit given that the most obscured AGN are missed in X-rays. We estimate the 
contribution of AGN missed in X-rays, using a population synthesis model, to 
be ~45% of the observed AGN contribution, making the AGN contribution to 
the infrared background at most ~2-10% in the 3-24 micron range, depending 
on wavelength, lower than most previous estimates. The AGN contribution to 
the infrared background remains roughly constant with source flux in the IRAG 
bands but decreases with decreasing flux in the MIPS 24 /rm band, where the 
galaxy population becomes more important. 

Subject headings: galaxies: active, infrared: galaxies 


1. Introduction 

Hard X-ray surveys are not strongly affected by dust obscuration and thus provide a 
relatively unbiased view of the AGN phenomenon. Deep Ghandra X-ray surveys have allowed 
us to test the basic prediction of population synthesis models of the X-ray background, 
namely, that a combination of obscured and unobscured AGN is needed to explain the X-ray 
background hard spectral shape (Setti & Woltjer 1989; Gilli et ah 2001). According to the 
AGN unihcation paradigm, obscuration comes from optically thick dust blocking the central 
engine along some lines of sight. The temperature in this structure, which can range up to 
1000 K (dust sublimation temperature), and the roughly isotropic emission toward longer 
wavelengths should make both obscured and unobscured AGN very bright in the mid-far 
infrared bands. 

This prediction was tested by IRAS (e.g.. Rush et ah 1993) and ISO (e.g., Fadda et al. 
2002; Matute et al. 2002; Franceschini et al. 2002) observations, which discovered a popu¬ 
lation of very luminous infrared sources coincident with X-ray emitters and therefore very 
good candidates to be obscured AGN. With relatively bright flux limits, however, these sur¬ 
veys sample the nearby Universe or extremely luminous, rare AGN. Very deep surveys were 
possible with ISO (e.g., Fadda et al. 2002), but only over a very small area. With the Spitzer 
Space Telescope, it is now possible to perform very deep infrared surveys, ~100-1000 times 
deeper than typical IRAS observations, over appreciable areas, ~10-1,000 times more area 
than the deep ISO surveys, thus permitting us, for the first time, to study a representative 
sample of the AGN population at cosmologically-signihcant distances. 

Since AGN comprise most of the X-ray background and are also very bright infrared 
sources, it is natural to ask what fraction of the extragalactic infrared background they 



3 


may contribute. Based on previous observations, the predicted contribution of AGN to the 
infrared background ranged from ~20% (Fabian & Iwasawa 1999) to just a few percent (Xu 
et ah 2001), depending mostly on the assumptions for the amount of energy re-radiated 
in the infrared and the spatial density of obscured AGN. The deeper infrared observations 
possible with Spitzer are important for understanding the true contribution of AGN to the 
extragalactic infrared background light. 

Multiwavelength observations of a large, uniformly selected sample of moderate-luminosity 
AGN at cosmologically significant distances provide a complete picture of the AGN emission 
across the frequency spectrum and establish the total contribution of these sources to the 
energy budget of the Universe. This is one of the main goals of the Great Observatories 
Origins Deep Survey (GOODS), a 0.1-square-degree study divided in two helds that overlap 
with the deepest X-ray observations, the Ghandra Deep Fields North (Brandt et al. 2001) 
and South (Giacconi et al. 2001). It consists of deep imaging in the optical with the Hubble 
Space Telescope (Giavalisco et al. 2004) and in the infrared with the Spitzer Space Tele¬ 
scope (Dickinson & Giavalisco 2002). These very deep multiwavelength observations allow a 
highly complete, relatively unbiased, view of the AGN phenomenon and a test of the AGN 
unihcation paradigm, well established in the local Universe, at higher redshifts. 

In this paper we present Spitzer observations of the AGN in the GOODS helds, from 
3.6 to 24 microns. In § 2 we describe the Spitzer IRAG and MIPS observations. Results 
and discussion are presented in § 3, addressing the infrared luminosities of AGN and normal 
galaxies, the X-ray to infrared ratios for GOODS and local AGN, the infrared number 
counts of AGN, and the AGN contribution to the extragalactic infrared background light. 
Our conclusions are presented in § 4. Throughout this paper we assume Hq = 70 km s 
Mpc“^, VLm = 0.3 and Da = 0.7. 


2. Spitzer Data and AGN Sample 

The GOODS North and South helds, which cover ~160 square arcminutes each, were 
observed with the Spitzer Space Telescope as part of the GOODS Legacy program (PI; 
Mark Dickinson). Observations were performed using the IRAG camera (Fazio et al. 2004b) 
in the four available bands (3.6-8/im) and the MIPS instrument (Rieke et al. 2004) in the 
24 /im band. Exposure times per sky-pixel for each held were ~23 hours in the IRAG 
bands. The MIPS 24-micron observations, available to us only in the GOODS-North held, 
had an exposure time of ~10.3 hours. Individual images were processed through the Spitzer 
Science Genter basic calibrated data pipeline, and then post-processed and combined using 



customized routines developed by the GOODS team^ (M. Dickinson et ah, in prep). Sources 
were detected using SExtractor (Bertin & Arnouts 1996) on a combined 3.6 and 4.5 /rm band 
image filtered using a 5x5 pixel Gaussian profile with a FWHM of 1.8". Photometry was 
performed using SExtractor. The AUTO magnitnde was used since aperture corrections, 
which we do not apply, would be very small, according to simulations^ (M. Dickinson et 
ah, in prep.). If instead we perform aperture photometry (with an apertnre diameter of 
4") and make the corresponding apertnre corrections, which range from 0.3 magnitudes in 
the 3.6 /im band to 0.7 at 8 /rm then the AUTO and apertnre-corrected magnitudes differ 
by less than 0.2 magnitndes on average, with no obvious dependence on magnitnde and no 
signihcant ontliers. Therefore onr results do not depend signihcantly on the method used to 
perform the photometry. 

The formal 5a, background-limited flux limits for an isolated point sonrce in the four 
IRAG channels are 0.13, 0.22, 1.45, and 1.61 /iJy respectively. In practice, the photometric 
uncertainty for a given source depends on the degree of crowding in the images, and on the 
details of how the photometry and sonrce deblending are carried ont. This increases the 
photometric errors, which we evaluate using Monte Garlo simnlations of artihcial sources 
added to the real data. According to simnlations, all the IRAG images are 90% complete at 
~14 /iJy. 

The MIPS 24-/rm image of the GOODS North held was rednced using a combination 
of the standard Spitzer Science Genter pipeline and customized routines created by the 
GOODS team (R. Ghary et ah, in prep.). Photometry was obtained at the IRAG-derived 
positions using a point spread function htting technique. The 5a hux limit for the MIPS 
observations is ~25 /iJy and are 84% complete at that hux. The relatively large aperture 
used to calculate huxes in both the IRAG and MIPS observations means that the light from 
the AGN and the host galaxy are snmmed; in any case, given the limited spatial resolntion, 
the two contributions cannot be separated from the infrared imaging alone. 

In order to dehne the AGN sample we start with the X-ray catalogs for the Ghandra 
Deep Fields North and South of Alexander et ah (2003), with the associations to the GOODS 
HST optical catalog from Bauer et ah (2004) which uses the likelihood method described by 
Bauer et al. (2000). Specihcally, our sample consists of sonrces detected in the hard (2-8 
keV) Ghandra band that do not have mnltiple optical counterparts in the GOODS AGS 


^More details about data reduction can be found at the webpage 
http://data.spitzer.caltech.edu/popular/goods/Documents/goodsjdataproducts.html 

^See the webpage http://www.noao.edu/noao/staff/ecm/simulations.html for more details and calcula¬ 
tion of aperture corrections and completeness simulations. 
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z-hand image. Only 3 sources in the South held and 6 in the North have multiple optical 
counterparts, so excluding these sources from the sample does not bias our results. (We do 
include sources with no z-hand counterpart down to the limit ;2~27 mag.) A large fraction 
of the sample, ~60%, have spectroscopic redshifts, obtained by Szokoly et ah (2004) in the 
South held and by Barger et ah (2003) in the North. Most of the X-ray sources are obscured 
AGN at moderate redshift, 0.5< 2 : <1, and unobscured AGN up to 2 : ~4. More distant 
obscured AGN are too faint for optical spectroscopy (Treister et al. 2004). Sources with 
spectroscopic redshifts and hard X-ray luminosity lower than 10"^^ ergs s“^ are probably not 
AGN dominated, since the most X-ray luminous star-forming galaxy known, NGG 3256, has 
a total X-ray luminosity Lx — 8 X 10"^^ ergs s ^ in the 0.5-10 keV band (Lira et al. 2002); 
so these were removed from our sample. This criterion removes 12 sources in the South held 
and 25 in the North. With these constraints, the hnal sample includes 126 X-ray sources in 
the GOODS-South held and 179 in the -North held. 

In order to look for Spitzer counterparts of the Ghandra sources, several search radii 
were tried, looking to maximize the number of sources with counterparts while keeping a 
low number of sources with multiple counterparts. The best results were obtained using a 
search radius of 2". Typical astrometric errors are ~0.5" in X-rays and ~0.2" in the IRAG 
bands; therefore, astrometric errors can be safely neglected when using a search radius of 
2". In the North held 160 sources (90%) have a single IRAG counterpart, while in the south 
held 118 sources (94%) have a single counterpart; the rest have multiple counterparts. This 
means that 100% of the X-ray sourees are also bright in the near-infrared. Most of the X- 
ray sources have an infrared counterpart separated by less than 1", and in many cases this 
distance is smaller than 0.5", as shown in Fig. 1. Assuming a Poissonian spatial distribution 
of infrared sources, the median chance of a random association is 0.74% and therefore we 
have high conhdence that the correct infrared counterpart was identihed correctly in all but 
one or two cases. The full catalog of AGN detected in the GOODS held used in this work is 
presented in the on-line version of the Journal, while for clarity a fraction of the information 
is presented in Table 1. 

Our sample only considers AGN selected based on the X-ray information. Because 
of this, highly obscured AGN in which most of the X-ray emission is absorbed will be 
missed in our sample. The ehects of this incompleteness are discussed in §4. Recently, 
AGN selection methods based only on IRAG and MIPS colors (e.g., Sajina et al. 2005) or 
a combination of Spitzer infrared and radio properties (e.g.. Mart Inez-Sansigre et al. 2005) 
were presented. In the work of Sajina et al. (2005), AGN are found in the region dehned by 
0.0Sog(^4.5/^3.6)<0.3, 0.2Sog(^5.8/^3.6)^.5, 0.3Sog(A8/^3.6)^.r, 0.lSog(A5,8/^4.5)^.3, 
0.3<log(S'8/>S'4,5)<0.5 and 0.1<log(S'8/*S'5,8)%0.3 using only the IRAG bands (Fig. 8 on Sajina 
et al. 2005). Although this criterion was dehned only for AGN at z<l, K corrections are 
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very small so this method can be extended to higher redshifts. Only 79 sources satisfy the 
AGN criterion of Sajina et ah (2005) in both the GOODS North and South helds. Of these, 
only 18 (~6%) sources are detected in hard X-rays. This means that using this selection 
method misses a large fraction of the AGN in the held, but it is useful to separate AGN 
continuum-dominated sources from normal galaxies. This can be explained because a power- 
law continuum is assumed for the infrared emission in AGN, while it is possible that most 
sources are dominated by the torus thermal emission in the infrared (Van Duyne et ah, in 
prep). 

The selection method of Martinez-Sansigre et al. (2005) is based mostly in radio emis¬ 
sion and it is effective to detect luminous obscured quasars. Using the infrared criterion of 
Martinez-Sansigre et al. (2005) on the GOODS North held, where MIPS data were available, 
only 3 sources have 5*24 >300 pJy and S'3.6<45 pJy and are detected in hard X-rays. Unfor¬ 
tunately, neither spectroscopic nor photometric redshifts are available for those sources, so 
we cannot (yet) classify them as obscured quasars. The fact that only a small number of 
sources are selected using this criterion is not surprising since only very luminous sources are 
found in this way and such sources are very rare and thus not easily detected in narrow-held 
surveys like GOODS. 


3. Results and Discussion 
3.1. Infrared Luminosity and Colors 

When spectroscopic redshifts were not available, we used photometric redshifts com¬ 
puted by Mobasher et al. (2004) for the South held and Barger et al. (2003) for the North 
held. Using both spectroscopic and photometric redshifts, we were able to calculate lumi¬ 
nosities for 93 sources (74%) in the South held and 105 sources (59%) in the North held. 
In order to study the quality of photometric redshifts we compared them with spectroscopic 
redshifts for sources in which both measurements are available. The average value of the 
quantity (2:spec-'2phot)/(^“*'^spec) is 0.02, which indicates that photometric redshifts can be 
used in a statistical sense to study the properties of a sample of sources. However, errors 
in the photometric redshifts for individual sources can be large. About 30% of the sources 
compared here had a diherence between photometric and spectroscopic redshift determina¬ 
tion greater than 0.3. Therefore, photometric redshifts are used only to calculate statistical 
properties of a sample, not for individual sources. 

The observed-frame luminosity distributions for GOODS AGN in the IRAG 3.6 pm and 
MIPS 24 pm bands are shown in Fig. 2. The vast majority of the sources can be classihed 
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as moderate-luminosity AGN, in the range 10^^ to ergs s“^. No K correction was 

attempted. However, from the theoretical AGN SED library of Treister et ah (2004), which 
spans unobscured to fully obscured AGN over 4 decades in luminosity, we estimate that 
infrared luminosities in the 3-20 micron range would change by at most a factor of two due 
to K corrections. This is smaller than the bin size used in these histograms, so we do not 
expect the hgure to change signihcantly were we able to do accurate K corrections. 

For comparison, in Fig. 2 we also show the observed-frame 3.6 and 24 micron luminosity 
distributions for normal galaxies in the GOODS fields, randomly selected to follow the same 
redshift distribution as the AGN sample. Redshifts for the normal galaxies were obtained 
from the spectroscopic surveys of Wirth et al. (2004) and Gowie et al. (2004) in the North 
held and from the VIMOS VLT Deep Survey (Le Fevre et al. 2004) in the South held. 
Gomparing the median infrared luminosities of AGN and normal galaxies, AGN are brighter 
by a factor of ~10 in the IRAG bands and ~15 in the MIPS 24 /rm band. There is a 
fair amount of overlap between the two populations, indicating that for low-luminosity AGN 
(~10^^-10^^ ergs s“^) the host galaxy can provide a signihcant contribution and perhaps even 
dominate the infrared emission. Additionally, the luminosity distribution for sources with 
hard X-ray luminosity lower than 10^^ ergs s“^ (and thus removed from our AGN sample) is 
also shown in Fig. 2. These sources have a luminosity distribution similar to that of normal 
galaxies and therefore we conclude that the X-ray emission in those sources is not from AGN 
origin and most likely from star-formation activity. 

In order to understand the relative contributions of infrared emission from active nucleus 
and galaxy, we performed a simple simulation in which we randomly paired a normal galaxy 
from the redshift-equivalent distribution (infrared emission from stars only) to each AGN 
in the sample (infrared emission from stars plus active nucleus). We then calculated the 
distribution of the ratio of galaxy to nuclear light, assuming the assigned galaxy correctly 
represented the emission from the host galaxy. This process was repeated 150 times with 
different random pairings. If indeed these normal galaxies correctly represent the host galaxy 
emission, which is of course uncertain but which is certainly true at optical wavelengths 
(e.g., Urry et al. 1999; Dunlop et al. 2003), the ensuing distribution indicates the relative 
importance of the two components. The results of a set of these simulations are shown 
in Fig. 3 for the 8 and 24 micron bands. We hnd that at 8 /xm, ~ 87% of the sample is 
dominated by the AGN (i.e., the luminosity of the active nucleus is greater than the galaxy 
luminosity). Separating the sample in three luminosity bins, each with equal numbers of 
AGN, we hnd that for total AGN luminosities below 4 x 10^^ ergs s“^, roughly 3/4 of the 
sample is AGN-dominated. For 4 x 10^^ < Lj^ < 1.6 x 10"^"^ ergs s“\ the fraction of AGN- 
dominated sources is roughly 90%, and for higher luminosities it rises to 97%. At shorter 
wavelengths and at 24 fim, the fractions are similar. At 24 /rm, ~80% of the sources are 



AGN-dominated in the lowest luminosity bin and ~90% at the highest luminosities. Of 
course, because of the limited spatial resolution, we cannot say from the infrared data alone 
whether the AGN host galaxies are similar to normal galaxies in the same redshift range. 
Given the similarity of the optical properties, however, it would require a fairly contrived 
scenario to have the infrared properties be massively mismatched. The correlation of infrared 
and X-ray flux (next Section) also supports the idea that the infrared emission from AGN 
is dominated by the active nucleus. 


3.2. X-ray to Infrared Ratios 

In Fig. 4 we show the ratio of hard X-ray (2-8 keV) to infrared (8 fim) luminosity versus 
hard X-ray and infrared luminosity for GOODS AGN with spectroscopic redshifts and thus 
reliable luminosities. Using the X-ray spectrum to separate obscured and unobscured AGN 
(the separation was set at Nh=10‘^‘^ cm“^, with Nh inferred from the X-ray hardness ratio, 
assuming an intrinsic photon slope of r=1.7), we hud no differences between the two classes 
of AGN (hlled and open triangles in Fig. 4). This is consistent with the AGN unihcation 
paradigm since both hard X-ray and infrared emission are relatively isotropic and unaffected 
by the amount of dust in the line of sight. Similar results are found if the 24 micron fluxes 
are used instead; however, these fluxes are only available for half of our sample and thus we 
decided to use the 8 micron fluxes. 

In order to compare with previously known sources, in Fig. 4 we also show the average 
observed-frame hard X-ray to infrared luminosity ratio for both local Seyfert 1 and 2 galaxies 
in three X-ray luminosity bins, from the compilation of Silva et ah (2004). The average ratio 
for the local sources is 0.10±0.02 (standard deviation for the sample, including K correction), 
considering both obscured and unobscured AGN and excluding the top and bottom 10% 
values of the distribution. For the GOODS sample this ratio is 0.18±0.12, also excluding 
the top and bottom 10% values, and there is also no signihcant difference in the average 
values for the obscured and unobscured AGN. The average value for the GOODS sample in 
three X-ray luminosity bins is shown by the open squares in Fig. 4. There is no statistically 
signihcant discrepancy between the GOODS and local samples, although there is a hint that 
the values might diverge at high X-ray luminosity (<2a signihcance). However, it should be 
noted that Malmquist bias, which would make the average luminosity to appear to be higher 
at higher redshifts, can play a role in this difference. 

The trend for the average ratio of X-ray to infrared luminosity to increase with X-ray 
luminosity in the GOODS sample can be explained by the contribution of the host galaxy 
to the infrared light at low luminosities. Assuming that this ratio has a value of 0.35 for the 
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highest luminosity sources, which are clearly AGN-dominated, and adding to the infrared 
luminosity the contribution from an average host galaxy (estimated from the median 8 /rm 
luminosity for our random sample of galaxies, 7.67x 10*^^ ergs s“^), we were able to reproduce 
the observed trend of X-ray to infrared luminosity ratio with X-ray luminosity (solid line in 
Fig. 4). 

The spread in the values of the X-ray to infrared luminosity for the GOODS sample can 
be explained by the combined effects of ignoring the K correction and the intrinsic spread in 
the luminosity of the host galaxy. Using the AGN spectra library of Treister et al. (2004), we 
estimate that the K correction can account for variations of ~ 50% around the mean observed 
ratio. The dashed lines in Fig. 4 show the effects of a spread in host galaxy luminosity, which 
is particularly important for low-luminosity AGN and can account for most of the observed 
spread in the X-ray to infrared luminosity ratio at these luminosities. 

Performing a partial correlation analysis we found no signihcant correlation between the 
hard X-ray to infrared luminosity ratio and either hard X-ray or infrared luminosity, or with 
redshift or absorbing column density {Nh)- We found only a mildly signihcant correlation 
between hard X-ray and infrared luminosities, with a Pearson coefficient of ~0.604 and a 
correlation of similar signihcance between hard X-ray and infrared huxes, indicating that 
the huxes in these bands are directly related to the luminosity of the central engine. The 
lack of correlation of the ratio with redshift or luminosity suggests there is little dependence 
of the emission mechanism on either redshift or luminosity. 


3.3. Infrared Counts of AGN 

The cumulative distributions of infrared huxes in the 3.6, 5.7, 8 and 24 micron bands for 
the GOODS AGN are presented in Fig. 5. The observed counts were corrected to account 
for the diherences in covered area depending on the X-ray hux (since the sample is X- 
ray selected) based on the relation given by Fig. 5 in Treister et al. (2004), obtained by 
adapting the sensitivity versus area function for the Ghandra Deep Fields of Alexander et al. 
(2003) to the portion covered by GOODS (D. Alexander and F. Bauer, priv. comm.). No 
correction for incompleteness in the X-ray or infrared catalogues is necessary as these are 
almost complete in the hux ranges sampled in the number counts determination. Error bars 
in these distributions were calculated following the prescription of Gehrels (1986) for 84% 
conhdence levels on the number of sources. The number counts in the 3.6, 4.5 and 5.7 yim 
bands can be well described by a power-law with integral slope 1.5±0.2 for huxes higher than 
50 /rJy from a maximum likelihood analysis of the unbinned data (Grawford et al. 1970). In 
the 8 fim band the integral slope is slightly shallower, 1.2±0.2, at huxes higher than 70 /rJy, 
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with a density of ~1000 AGN per square degree at a depth of ~35 /xJy. In the MIPS 24 /rni 
band the slope is shallower still, 0.9±0.2, at all the fluxes sampled by the GOODS sources, 
and there are 1000 AGN per square degree at a flux density of ~200 /xJy. The resulting best 
hts are shown in Fig. 5. In order to test the goodness of these hts we performed a K-S test 
comparing the observed distribution to the fitting function, finding that the null hypothesis 
(that both distributions are drawn from the same parent distribution) cannot be rejected 
at conhdence levels higher than ~40% in all the bands. The values of the K-S probability 
for each band are 45%, 50%, 75%, 85% and 40% on the 3.6, 4.5, 5.7, 8 and 24 /rm bands 
respectively. 

In the unihcation paradigm for AGN, most of the energy absorbed by the circumnuclear 
optically thick matter is re-emitted at infrared wavelengths. It is not surprising then that 
AGN are luminous infrared sources, even when the optical, ultraviolet and soft X-ray emission 
is faint. Population synthesis models for the X-ray background suggest that a substantial 
fraction of AGN are obscured along the line of sight, thus will be bright only at mid- to 
far- infrared and hard X-ray wavelengths. For a simple unihcation model, in which the 
material is distributed in a torus and the ratio of obscured to unobscured AGN is set to 
~3:1, as observed in the local Universe (Risaliti et ah 1999), we can calculate the expected 
infrared counts. Model details are as described in Treister et al. (2004) and modihed by 
Treister & Urry (2005) to include a decrease in the number of obscured AGN with increasing 
luminosity, with the exception that we no longer interpolate between the near-infrared and 
optical spectra because this artihcially over-predicts the mid-infrared hux. 

The obscured-to-unobscured AGN ratio for sources with spectroscopic redshifts (and 
thus Nh determinations) in our sample is ~1.5:1, but this is consistent with the assumed 
ratio in the model, given that most of the obscured AGN do not have spectroscopic redshift 
since they are optically faint, and some obscured AGN are even missed in X-rays. Therefore, 
the observed value of the ratio only represents a lower limit to the intrinsic ratio, which can 
easily be higher by a factor of two and therefore consistent with the average value assumed 
in our model. 

This simple AGN unihcation model agrees within a factor of two with the observed 
distribution in the GOODS helds over the IRAG 3.6-8 micron and MIPS 24 micron bands, 
as shown in Fig. 5. Gompared to the model, which assumes only an elliptical host galaxy, 
with a single (L*) luminosity and no star formation, and hxed parameters for the dust torus, 
there are slightly fewer observed AGN at the bright end (>100 /iJy). The small area of the 
GOODS survey causes bright AGN to be poorly sampled, hence the large error bars and the 
marginal signihcance of the discrepancy. Additionally, it is important to point out that in a 
cumulative plot the bins are not independent. Wide area surveys like SWIRE (Lonsdale et 
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al. 2003) and the MIPS GTO^ programs will provide a better estimate of the AGN counts 
at the bright end. At the faint end (<10 /xJy), the observed number of AGN is higher than 
the model predictions except in the 3.6 /xm band where there is very good agreement. These 
discrepancies are ~13% at 4 /rm and rising to 36% at 8 /xm and ~30% at 24 /rm. These 
discrepancies may suggest that for fainter AGN, the host galaxy plays a relatively more 
important role in some AGN, thus making them brighter. At 3.6 yum the stellar component 
in the host galaxy dominates; the good agreement with our simple model means that the 
host galaxy must be reasonably well represented by our assumed L* elliptical galaxy with 
no recent star formation. 


3.4. Extragalactic Infrared Background Light 

Using the Spitzer data we can estimate the integrated contribution of AGN to the 
extragalactic infrared background light at IRAG and MIPS wavelengths. A representative 
lower limit can be obtained simply by summing the infrared emission from the AGN detected 
in the GOODS helds, weighted by the effective area at that X-ray flux. In Fig. 6, we show the 
resulting contribution measured from the GOODS sources, compared to the observational 
limits for the total extragalactic infrared background compiled by Hauser & Dwek (2001; 
the lower limit is from the integrated extragalactic light, while the upper limit comes from 
fluctuations analysis) and updated measurements at the IRAG and MIPS 24 /rm wavelengths 
(Matsumoto et al. 2005; Fazio et al. 2004a; Kashlinsky et al. 1996; Stanev & Franceschini 
1998; Papovich et al. 2004). AGN and their host galaxies contribute approximately 
3.01^;3%, 5.6l3;9%, in the IRAG 3.6, 4.5, 5.7 and 8 jum bands respectively, and 

10.7;“,“% in the MIPS 24-micron band. (The uncertainties are dominated by the range in 
the total infrared background light.) These are lower limits to the total AGN contribution, 
since they exclude the contribution from faint sources below the GOODS X-ray detection 
limit (including a large expected population of Gompton-thick AGN; Treister & Urry 2005) 
and from bright sources that are not well sampled in the relatively small GOODS helds. 
Using our AGN unihcation model, we can calculate the relative contribution of these X- 
ray faint sources and add it to the total emission from X-ray detected sources, in order to 
estimate the total contribution of AGN to the infrared background. We hnd that sources 
missed in X-rays can add ~45% more hux in the 3-24 micron range, with ~30% extra hux 
coming from obscured AGN that are too faint to be detected in X-rays and the remaining 
~15% coming from bright sources that are not well sampled in the GOODS helds. Gorrecting 
for this factor, the AGN contribution is increased to ~3-8% in the IRAG bands and ~15% in 


^More information can be found at the web site http://lully.as.arizona.edu/ 
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the MIPS 24-micron band. Given that these values are obtained from the summed emission 
of AGN and their host galaxies, these are upper limits to the AGN contribution to infrared 
background, or equivalently, can be seen as the contribution to the infrared background by 
AGN plus their host galaxies. 

We also investigate the flux dependence of the AGN contribution to the integrated 
infrared emission. This was computed by comparing the total flux from AGN with the 
integrated flux from all sources for several flux bins. At 3.6 microns the AGN contribution 
is 2.1% above 600 juJy, rising to ~4% above 10 juJy. At 8 microns, the AGN contribution 
is ~6% and remains mostly constant, with a slight increase to ~10% at 40 pJy. Only in 
the MIPS 24 pm band do we see a decrease in the AGN contribution with decreasing flux, 
going from 30% above 1000 pJy to ~8% above 100 pJy, as shown in Fig. 7. In other words, 
at lower fluxes the number of galaxies rises much more steeply than the number of AGN. 
While this can also mean that many faint AGN are missed by the X-ray selection, in Fig. 7 
we use the previously described AGN model to estimate the AGN fraction as a function of 
infrared flux, correcting for AGN missed by the X-ray selection. As can be seen in Fig. 7, 
the decrease in the AGN fraction with decreasing infrared flux is still observed, even after 
correcting for AGN missed in X-rays. 

Previous reports of the AGN contribution in general estimated a higher contribution to 
the integrated emission. Earlier infrared observations with ISO (e.g., Fadda et ah 2002) in 
the 15-micron band at brighter fluxes give higher values for the AGN contribution (17% at 
15pm) but depend on large extrapolations to lower fluxes. (The GOODS MIPS observations 
are ~25 times deeper than the ISO data on the Lockman Hole.) Matute et ah (2002) made 
an estimate based on ISO observations of bright unobscured AGN in the ELAIS SI field. If, 
in order to calculate the total AGN contribution from their measured value for unobscured 
AGN we assume an obscured-to-unobscured ratio of ~4:1, their results are consistent with 
the observed values obtained from the GOODS observations. An estimate based on X-ray 
samples and an average X-ray-to-infrared ratio by Silva et al. (2004) is ~50-80% higher 
than the values obtained directly from the Spitzer observations. Using a backward evolution 
model based on IRAS and ISO observations and infrared luminosity functions, the AGN 
contribution to the infrared background was also calculated by Xu et al. (2001). In this case 
only an infrared color criterion was used to identify AGN. This prediction is similar to the 
observed values in the 3.6 and 4.5 pm bands but estimates a lower contribution at longer 
wavelengths, ~55% lower at 24 pm. A recent study of AGN detected by Ghandra using 
Spitzer performed by Franceschini et al. (2005), reported a value for the AGN contribution 
to the infrared background of 10%-15% at 24 pm in good agreement with our results. 


The earlier estimates of the AGN contribution to the integrated extragalactic infrared 
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light are therefore largely consistent with our results; the earlier numbers differ mainly 
because of large extrapolations to lower fluxes. Summarizing, the more precise values from 
our deep Chandra and Spitzer observations are 2.4% at 3.6 pm, 3.0% at 4.5 pm, 5.6% at 
5.7 pm, 4.7% at 8 pm and 10.7% at 24 pm. These numbers should be multiplied by a factor 
of ~1.45 in order to account for the contribution by Compton-thick AGN not detected in 
X-rays by Chandra, namely, 3.5% at 3.6 pm, 4.4% at 4.5 pm, 8.1% at 5.7 pm, 6.8% at 8 pm 
and 15.5% at 24 pm. 


4. Conclusions 

We present here Spitzer observations of the moderate-luminosity, moderate- to high- 
redshift AGN in the GOODS helds, which include the deepest infrared and X-ray obser¬ 
vations to date. In general these sources have similar X-ray to infrared luminosity ratios 
as AGN in the local Universe. AGN are ~10-15 times brighter in the infrared than normal 
galaxies, although there is some overlap between the luminosity distribution of AGN and nor¬ 
mal galaxies, indicating that for low-luminosity AGN the contribution to the infrared light 
from the host galaxy can be very important and may well dominate the observed emission. 
This can explain the observed trend of the X-ray-to-infrared luminosity ratio to increase 
with X-ray luminosity. 

We report the observed infrared flux distribution for GOODS AGN, which has a slope 
of ~1.5 to 0.9 for fluxes higher than ~100 pJy, with a source density of ~1000 AGN per 
square degree at a depth of ~50 pJy in the IRAC bands. Comparing to the expectations of a 
unihcation model for AGN, we hnd reasonable agreement, showing that the basic properties 
of the moderate to high redshift AGN population are well explained by the same unihca¬ 
tion paradigm that hts local sources very well, but with a factor of two less dust in the 
circumnuclear region. 

From the Spitzer observations, we calculate the minimum AGN contribution to the ex- 
tragalactic near-mid infrared background, obtaining a lower value than previously estimated, 
ranging from 2% to 10% in the 3-24 micron range. Accounting for heavily obscured AGN 
that are not detected in X-ray, the AGN contribution to the infrared background increases 
by ~45%, to ~3-15%. 

We thank the GOODS team for very useful comments and the anonymous referee for 
suggestions that improved the paper. ET thanks the support of Fundacion Andes, Centro 
de Astrofisica FONDAP and the Sigma-Xi foundation through a Grant in-aid of Research. 
This work was supported in part by NASA grant HST-GO-09425.13-A. 
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Table 1. Catalog of AGN in the GOODS fields. 


Field 

ID 

Prob.*^ 

X-ray Flux (erg 
Soft (0.5-2 keV) 

cm ^ s 

Hard (2-8 keV) 

Hardness Ratio^ 

3.6 /^m 

Infrared Flux (^Jy) 

4.5 fim 5.7 fim 8 fim 

24 fim 

Redshift 

Redshift type^ Classification*^ 

S 

34 

0.011 

3.58e-15 

7.03e-15 

-0.42 

24.89 

26.55 

30.76 

33.73 

-1.00 

1.037 

s 

B 

S 

39 

0.007 

2.77e-15 

5.62e-15 

-0.41 

8.24 

8.47 

9.55 

14.32 

-1.00 

2.726 

s 

B 

s 

41 

0.004 

4.35e-15 

8.61e-15 

-0.42 

52.97 

60.26 

59.70 

53.95 

-1.00 

1.616 

s 

B 

s 

42 

0.032 

2.04e-16 

6.74e-16 

-0.24 

5.06 

6.49 

9.46 

5.81 

-1.00 

9.999 

N 

U 

s 

43 

0.015 

6.84e-16 

8.99e-16 

-0.54 

7.18 

7.94 

6.25 

6.73 

-1.00 

2.573 

s 

N 


^'Probability of random association assuming a Poissonian distribution of sources. Expression l-exp(-n tt d?) was used where n is the density of sources and d is the distance 
between the X-ray and infrared sources. 

^Defined as (H-S)/(H-hS) where H and S are the background-substracted counts in the Hard and Soft X-ray bands respectively. 

^S: Spectroscopic Redshift, P: Photometric Redshift, N: No redshift measured 
Broad line AGN, N; Narrow line AGN, G; Galaxy, U: Unknown 


Note. — This table is published in its entirety in the electronic edition of the Astrophysical Journal. A portion is shown here for guidance regarding its form and content. 








18 



X-ray/IR Offset (arcsec) 


Fig. 1.— Distribution of offsets between the X-ray and infrared positions for sources in our 
AGN sample. Most of the X-ray sources have a infrared counterpart closer than 1", while 
for a large number of AGN, the distance between the X-ray and infrared positions is smaller 
than 0.5". Therefore, we are very conhdent that the right infrared counterpart for the X-ray 
sources in our sample was used in all but one or two cases. 
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Fig. 2.— Infrared luminosity distribution for AGN with spectroscopic {dashed lines) and both 
spectroscopic or photometric redshifts {solid lines) in the IRAC 3.6 pm {left panel) and MIPS 
24 pm {right panel) bands. Most of the sources in our sample can be classihed as moderate- 
luminosity AGN, in the ergs s“^ range. Dotted lines: Luminosity distribution for 

random GOODS held galaxies with either photometric or spectroscopic redshifts selected 
to have the same redshift distribution as the AGN sample. Sources classihed as AGN in 
X-rays are ~10-15 times brighter in the infrared than normal galaxies. Dot-dashed/hatched 
histograms: Luminosity distribution for X-ray sources with Lx <10"^^ ergs s“^, and thus not 
included in our AGN sample. These sources have distributions similar to normal galaxies 
and much fainter on average than the AGN luminosity distribution. Therefore, we conclude 
that the X-ray emission in these sources does not come from AGN. 
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Fig. 3.— Simulated fraction of AGN to host galaxy infrared luminosity as a function of total 
infrared luminosity in the 8 /rm [left panel) and 24 gm [right panel) bands. The AGN and 
host galaxy emission are separated by associating the luminosity of random galaxies in the 
held to AGN in our sample. Here the results for 5 random associations for each source are 
shown. As can be seen, most of the sources are AGN dominated in the infrared with the 
fraction increasing with increasing total luminosity. This result does not depend much on 
the observed wavelength. 
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Fig. 4.— Left panel: Hard X-ray (2-8 keV) to infrared (8/im) luminosity ratio as a function of 
hard X-ray luminosity for obscured {solid triangles) and unobscured {open triangles) AGN 
(using Nh=10‘^^ cm“^ as the separating point) with measured spectroscopic redshifts in 
the GOODS fields. Average values (excluding the top and bottom 10%) in three X-ray 
luminosity bins are shown by the open squares, with error bars corresponding to the RMS 
dispersion. For comparison, the average values for this ratio obtained from a compilation of 
local AGN by Silva et al. (2004) are also shown {filled eircles: obscured AGN, open eireles: 
unobscured AGN). GOODS AGN, the majority of them at ~0.5-l (higher-redshift AGN 
are too faint for spectroscopy), and local sources have similar values for this ratio, possibly 
diverging at higher luminosity. Solid line: X-ray to infrared luminosity ratio assuming an 
intrinsic ratio of 0.35 and a host galaxy with a luminosity of 7.67x10^^ ergs s“^ (median 
luminosity for our random sample of GOODS normal galaxies). The dashed lines show the 
effects of changing the luminosity of the host galaxy by one order of magnitude (lower line: 
brighter host galaxy). The small observed trend with luminosity can be explained by the 
relatively more important contribution from the host galaxy for low-luminosity AGN. Right 
panel: Hard X-ray to infrared luminosity ratio as a function of infrared luminosity. Symbols 
are the same as in the left panel. In this case, sources below the X-ray luminosity threshold 
are also included in the diagram. The same conclusions are still valid when the X-ray to 
infrared luminosity ratio is measured as a function of infrared or hard X-ray luminosity. 
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0.01 0.1 1 0.01 0.1 1 

Flux (mJy) 

Fig. 5.— Filled circles: Observed cumulative infrared flux distributions for X-ray detected 
AGN in the GOODS North and South fields in three IRAG bands (3.6-8/im) and for the 
North held only in the MIPS 24yum band. Error bars correspond to the 84% conhdence level 
in the number of sources on that bin (Gehrels 1986). Note that because this is a cumulative 
plot, the bins are not independent. Solid line: Power-law ht to the observed number counts. 
Dashed line: Expected infrared hux distribution based on the models of Treister et ah (2004), 
taking into account the X-ray hux limit in the GOODS helds and correcting for the ehects 
of undersampling of the bright end caused by the small volume of the GOODS helds. Dotted 
line: Expected infrared hux distribution not considering the X-ray hux limit (i.e., all sources 
expected from extrapolating the AGN luminosity function). In general, the model follows 
the normalization and shape of the counts reasonable well, especially considering the poor 
statistics at the bright end. At the faint end, where the host galaxy light dominates, the 
model diverges more siernihcantlv since the model assumntions are clearlv far too simnle. 
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Fig. 6.— Extragalactic infrared background intensity as a function of wavelength. Shaded 
region-. Allowed values compiled by Hauser & Dwek (2001). 1-4 /rm: Measurements reported 
by Matsumoto et al. (2005). 4-8/rm: Lower limits from Fazio et al. (2004a), upper limits 
from Kashlinsky et al. (1996) at 4.5 /rm and Stanev & Franceschini (1998) at 5.8 and 8 /rm. 
Measurement at 24 /rm from Papovich et al. (2004). Solid eireles: Integrated infrared emis¬ 
sion from X-ray detected AGN in the GOODS helds. Error bars in these measurements are 
estimated based on the number of sources. Triangles: Integrated AGN emission from the 
models of Treister et al. (2004) and Treister & Urry (2005) if only X-ray detected AGN are 
considered {lower) and including all AGN {upper). Dotted lines: Expected AGN integrated 
emission from the AGN models of Silva et al. (2004; upper line) and Xu et al. (2001; lower 
line). Open sguares: Expected AGN integrated emission at 15/im from the extrapolation to 
fainter fluxes of IRAS and ISO (Matute et al. 2002; M02; corrected by a factor of 5 to include 
the contribution from obscured AGN) and ISO only (Fadda et al. 2002; F02) observations. 
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Fig. 7.— Fraction of sources classified as AGN increases with infrared flux in the 24 /iui 
band {solid circles). Horizontal error bars show the size of the flux bins, while vertical error 
bars show the la Poissonian errors on the number of sources. In order to test if this effect 
can be caused by the X-ray selection, we also show the contribution corrected by the AGN 
expected to be missed by X-ray selection, as estimated using the AGN population synthesis 
model {open circles). As can be seen, this effect is still observed after correcting for the AGN 
not detected in X-rays. 



